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Introduction 


■ Lattice Boltzmann (LB) Method is a relatively new method for 
flow simulations. 

■ The start point of LB method is statistic mechanics and 
Boltzmann equation. 


■ The LB method tries to set up its model at molecular scale 
and simulate the flow at macroscopic scale 

■ LBM has been applied to mostly incompressible flows and 
simple geometry 
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Statistic Mechanics 




Statistic mechanics views fluid as 
a collection of particles. 

The properties of the fluid are 
determined by the average 
properties of the particles in the 
collection 


A small domain near point x 
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Boltzmann Equation 

1) Distribution Function 



Phase space for ID problem 
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A phase space consists of 
both location and velocity of 
particles is introduced. 

The distribution function is 
defined as the density of the 
number of particles at point 
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Boltzmann Equation 

2) Time Evolution Of Distribution Function 


at time t + At 









Assume that there is no 
external force and no 
collision between the 
particles, the velocity of the 
particles will not change 


The two domain have 
the same number of 
particles 
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f(x + £At,£,t + A t)dxd£ = f(x,£,t)dxd£ 

f(x + %A t, g, t + At) = f(x, g, t) 
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Boltzmann Equation 

3) Collision Term 


at time t + At 



— >■ 
x 


Collisions between particles 
change their velocities, and make 
them move in and out of the 
domain 

A collision term describes the 
net increase of the density of the 
number of particles in the 
domain due to the collision 



COMPUTING 
INFORMATION & 
COMMUNICATIONS 
TECHNOLOGY 
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Boltzmann Equation 

4) BGK Collision Term 


■ One of the simplest collision models is the Bhatnagar, Gross and 
Krook (BGK) simplified collision model 


T 

■ The BGK model is widely used in LB models 
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Boltzmann Equation 

5) Boltzmann Equation 


The Boltzmann equation in finite difference form: 

/(x + %Ai,%,t + Af) = £2At 

By Taylor expansion, the above equation can be written in the 
differential form 

f -+<f-v/ = n 

Ot 


This is the Boltzmann Equation 
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The Macroscopic Properties 


Y (x t)- ^ 

i (!)/(*> 

I yX 5 f ) 



The Macroscopic properties 
are determined by the 
average value of properties 
of the particles 


density 

i o = mjf(x,^,t)d 3 ^ 

momentum 

pv ; = mjgf(x,g,t)d 3 g 

thermal energy 
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The Enskog-Chapman Expansion 




It has been shown that the Euler equation and Navier-Stokes 
equation are the zeroth-order and first order approximations 
of the Boltzmann equation, respectively., 

That is to say, the Boltzmann equation describes the fluid 
phenomena in a more accurate way than tradition fluid 
dynamics does 
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Lattice Boltzmann Method 


■ Lattice Boltzmann Method can be reviewed as a numerical 
method to solve the Boltzmann equation 

■ In LB method, the phase space is discretized. 

■ In a LB model, the velocity of a particle can only be chosen 
from a velocity set, which has only a finite number of velocities 
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Lattice Boltzmann Method 

1) Velocity Set and Grid 


1 
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For the convenience of 
computation, the position 
space is discretized in such a 
way that the particles travel 
with one of the velocity in the 
velocity set will arrive at a 
correspondent node at next 
time step. 
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The Lattice Boltzmann Method 

2) Macroscopic Properties 


■ In the LB method, the macroscopic properties are evaluated 
through the weight summation 


y = j 7i{z)f M [x,zM = I 

a 



X 

ii 

M 

if"* 

S-x 

171 1 p 

P£ = ^2ji*- v 

2 f{eq) 
J a 

a 

a 

2 „ 1 
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f ( : q) i 


\=WJ M \ 
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Lattice Boltzmann Method 

3) Time Evolution Equation 


Boltzmann equation in finite difference form with the BGK collision 
term 


/(* + la J + A? ) = /(*>£, . <) (/fa, la . t)~ f ieq) (*>!«><)) 

T 

where P is one of velocities in the velocity set 

z>a 


Set A t = T 


/I 

[x + £a At ’£a’ t + Af) 

II 

(-*■’ ^>a ’ 
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The Lattice Boltzmann Method 

4) The Calculation 


a 

[x-la A t,l a ,t-At) 

i=2Xi 

a 

[x -£ a At,t - At) 


■ The calculation of LB method can be viewed as following steps 

1) a node dividing its macroscopic quantities into parts 

2) the node sending parts of the macroscopic quantities to corresponding 
nodes 


3) a node receiving the parts of the macroscopic quantities sent by 
nearby nodes, adding them together and obtaining the macroscopic 
quantities for next time step 
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Advantages of LB method 


■ Start from a clear and direct physical picture at molecular 
level; 

■ Algorithm is simple and straightforward; 

■ Natural parallel scheme; 






■ Easy to incorporate the physical phenomena at molecular 
level, possible of modeling fluid phenomena that can not be 
modeled with the traditional fluid mechanics. 
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Motivation 


■ There are very few reports about the successful LB 
simulations of the real life problems 


■ Successful extension of LB method to the simulation of 
turbomachinery can show the maturity of LB method and the 
promise of this method in the simulation of real life problem 






■ The parallel nature of LB method make it a possible high 
performance solver for turbomachine simulations 
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Challenges 


1. Most of LB models can only simulate the flow with a small 
Mach number 

2. Most of past LB simulations are laboratory type simulations 
with simple computational domain and boundary 


3. The proper LB model for simulation should be developed 
and necessary techniques should be developed 
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The Current Work 


■ Develop a compressible LB model. 


■ Introduce a boundary condition that allows accurate turbine 
blade simulation 

■ Develop mesh treatment for the irregular computational 
domain 
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The Current Work 


■ Simulations carried out for three different cascades. It is the 
first time that turbomachines have been simulated by a LB 
model. 


■ The parallel performance of the LB model has been tested 
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Current Model 

1) The Velocity Set 
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■ A velocity in the velocity set 
consists of three parts. 

■ First the macroscopic velocity 
has been included explicitly into 
the microscopic velocity. 

■ Second, there is a set of 
diffusion velocity 
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Current Model 

2) The Diffusion Velocity 

■ There are 3 levels of diffusion 
velocity in current model 




c 


22 


— / 


c 


32 







C 


42 


V 


The module of the diffusion 
velocities are determined by 
macroscopic quantities 
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Current Model 

3) The Velocity Set 



of 


v,. 


IS 


A third se 
velocities 
introduced to carry 
the particles to the 
nearest vertex 
nodes 


r jvk = v + c j v + % 
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Current Model 

4) The Microscopic Quantities 


■ Mass, only one species is considered, is a constant 

m - 1 

■ Velocity, consists of macroscopic velocity and a diffusion velocity 

£jv — V Cjv 

■ Total energy is the same for all particles 

. 1 2 

c = — v +e 
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Current Model 

4) Modification of Microscopic Quantities 


■ For the purpose of recovering the correct Navier-Stokes equation, a 
correction term has been introduced 

p D 

Zjrk = 


p~b 0 d 0 2c 


7 T&v-v'k) 


Mass 

Velocity 

Total energy 


=1 “ V 


Zi* =v + c' v -/l M v 


jv jvk 


ijvk - ^jvk )f 2 


A 


v +e 
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Current Model 

5) The Equilibrium Distribution Function 

4 

M — rM 


f{eq) _ f 
J jvk j vk 

fi eq) =oc t d 


J _ (p-VoK 2 -D(y-\)pe 
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The Flow Chart Of Computation 
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The Mapping Between Coordinates and Indexes 




n 


(x,y)- 


(x,y) 


Mapping XY To MN 



■( m,n )■ 




Mapping MN To 
XY 


(m,n) 

The index of array 




The coordinate of nodes 
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Implement of Boundary Condition 

1) Wall Boundary - Auxiliary Nodes 




Auxiliary 


nodes 



Auxiliary nodes were 
introduced to implement wall 
boundary conditions 


Computational 

domain 
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Implement of Boundary Condition 

2) Wall Boundary - Extrapolation 
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Computationa 
1 Nodes 


Extrapolation 

Point 




Normal of 
Boundary 


The macroscopic 
properties of 
auxiliary nodes are 
extrapolated from 
the computational 
domain. 


Extrapolation 

Nodes 
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Implement of Boundary Condition 

3) Periodical Boundary 


Periodical Boundary 


Periodical Boundary 


Modified Destination 


i i i i i i i 


i i i i i i i 
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Wedge Cascade 

Geometry 





Wedge Cascade 

Result: Stream Line and Ma Contour 
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| C3X Cascade 

Geometry and Inlet and Outlet Parameters 







Stagger Angle 

59.89° 



Chord 

14.49cm 



Spacing 

11.77cm 



Solidity 

1.23 



Axial Chord 

7.82cm 


X 


Run Number 

a,n 

Pn (Pa) 

T n (K) 

Ma j 

Rej 

PilPa 

RUN 143 

0 

7755 

811 

0.17 

0.63 x 10 6 

0.50 

RUN 144 

0 

7889 

815 

0.16 

0.63 xlO 6 

0.59 
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C3X Cascade RUN 143 

esults 



C3X Cascade RUN 143 


Results 



P/Po 1 

1.00 

0.96 

0.92 

0.88 

0.84 

0.80 

0.76 

0.73 

0.69 

0.65 

0.61 

0.57 

0.53 

0.49 

0.45 



RUN 143 



2002 CISO Review 


NASA/TM— 2003-21 1896 


168 





i 


C3X Cascade RUN 143 

Results 


p/p 


ti 
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C3X Cascade RUN 144 

Results 


2002 CISO Review 




2002 CISO Review 


NASA/TM— 2003-21 1896 


169 


C3X Cascade RUN 144 
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C3X Cascade RUN 144 

Results 


p/p 


ti 
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VKI Cascade 

Geometry and Inlet and Outlet Parameters 



Stagger Angle 
Chord 
Spacing 
Solidity 
Axial Chord 


55.0° 

67.646 mm 
57.50 mm 
1.1765 
36.98 mm 


Run Number 

<Xtn 

Pa (Pa) 

r, (K> 

Ma j 

Re! 

Ma 2 

MUR 129 

0 

18200 

409.20 

0.15 

0.27 x 10 6 

0.84 
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VKI Cascade MUR 129 

Results 




2002 CISO Review 


NASA/TM— 2003-21 1896 


171 


VKI Cascade MUR 129 



± 


VKI Cascade MUR 129 

Results 


p/pii 
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Parallel Computing 


■ The LB method is a natural parallel method; 

■ The LB model is explicit in time; 

■ The area of dependent domain of a node is determined by the 
magnitude of velocity set, which is usually a small number. 
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Parallel Computing 

1) Division of Blocks 


■ The simulation of wedge cascade was parallelized to test the 
parallel performance of current LB model 
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Parallel Computing 

2) Information Exchange 
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Conclusion 


1) A compressible LB model has been successfully developed for 
turbomachinery simulations. 

2) Successful simulation of cascades has been carried out and it 
is the first successful turbomachinery simulation by a LB 
model. 


COMPUTING 
INFORMATION & 
COMMUNICATIONS 
TECHNOLOGY 


2002 C/SO Review 


4 


Conclusion 


3) A treatment of boundary condition in LB method has been 
introduced to the current compressible LB model. 

4) A new mesh treatment method has been devised in order to 
use regular mesh on a irregular geometry. 
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Conclusion 


5) The parallel efficiency of the new compressible LB model is 
studied. A linear efficiency has been demonstrated. 

6) The theoretical basis of the current model is analyzed in 
detailed. 
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